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CONVERSION  FACTORS.  NON- SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


_ Multiply _ 

degrees  (angle) 
inches 

miles  (US  statute) 
pounds  (force) 

pounds  (force)  per  square  inch 
pounds  (mass)  per  cubic  foot 
pounds  (mass)  per  cubic  inch 
square  inchc"; 


_ ^ _ 

0.01745329 

2.54 

1.609347 

4.448222 

6.894757 

16.01846 

27.6799 

6.4516 


_ To  Obtain _ 

radians 

centimetres 

kilometres 

newtons 

kilopascals 

kilograms  per  cubic  metre 
grams  per  cubic  centimetre 
square  centimetres 
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STRENGTH  PROPERTY  ESTIMATION  FOR  DRY.  COHESIONLESS 


SOILS  USING  THE  MILITARY  CONE  PENETROMETER 


PART  I:  INTRODUCTION 

Background 

1.  The  US  Army  Engineer  Waterways  Experiment  Station  (WES)  has  success¬ 
fully  used  cone  index  (Cl)  values  obtained  from  the  Military  Cone  Penetrometer 
Test  (MCPT)  (see  Figure  1)  to  characterize  soil  strength  since  the  1950's. 

Over  the  ensuing  AO  years,  WES  accumulated  a  large  database  of  Cl  values  from 
around  the  world.  These  Cl  values  are  used  in  empirical  vehicle-soil  interac¬ 
tion  relationships  to  predict  vehicle  performance  (Rula  and  Nutall  1971) .  The 
empirical  relationships  have  been  developed  and  refined  based  on  extensive 
field  testing  of  tracked  and  wheeled  vehicles.  However,  these  empirical  rela¬ 
tionships  are  limited  in  their  ability  to  model  the  detailed  physics  of 
vehicle-soil  interaction  required  for  describing  the  performance  of 


Figure  1.  Standard  military  cone  penetrometer 
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maneuvering  vehicles,  the  deformation  and  compaction  of  soils  resulting  from 
traffic,  and  next -generation  vehicles  whose  characteristics  fall  outside  the 
range  of  the  empirical  relations.  To  address  this  problem,  vehicle  -  terrain 
interaction  models  have  become  increasingly  theoretical,  relying  on  Newtonian 
physics  and  complicated  track-soil  and  tire-soil  interaction  models  to  predict 
vehicle  performance.  The  newer,  more  theoretical  models  require  that  the  soil 
be  characterized  by  fundamental  engineering  properties  such  as  cohesion,  fric¬ 
tion  angle,  shear  modulus,  and  unit  weight  rather  than  simple  index  properties 
such  as  Cl . 

2.  The  required  soil  properties  could  be  obtained  directly  from  a  tri- 
axial  shear  or  direct  shear  test;  however,  tests  such  as  those  are  expensive 
and  time-consuming  and  could  not  begin  to  replace  the  existing  body  of  accumu¬ 
lated  Cl  data.  The  MCPT,  on  the  other  hand,  is  a  simple,  easy-to-operate , 
man-portable  device  that  provides  a  rapid  and  inexpensive  means  of  indirectly 
gaging  soil  strength  in  the  field.  Unfortunately,  the  MCPT  only  provides  an 
indirect  measure  of  soil  strength  (the  index  property  C' )  rather  than  a  direct 
measure  of  strength  and  deformation.  In  order  to  make  use  of  both  the  exist¬ 
ing  and  future  bodies  of  Cl  data  worldwide,  correlations  with  which  we  car, 
infer  engineering  properties  from  Cl  measurements  must  be  developed. 

3.  A  first  attempt  at  inferring  engineering  properties  from  MCPT  data 
was  undertaken  by  Meier  and  Baladi  (1988).  Their  methodology  is  based  on  a 
theoretical  formulation  of  the  Cl  problem  using  cavity  expansion  theory  to 
relate  cone  penetration  resistance  to  the  strength  and  deformation  properties 
of  the  soil.  The  cavity  expansion  formulation,  originally  derived  by  Rohani 
and  Baladi  (1981),  incorporates  three  mechanical  properties  (cohesion,  fric¬ 
tion  angle,  and  shear  modulus)  and  the  total  unit  weight.  Obviously,  these 
four  unknown  soil  properties  cannot  be  back-calculated  directly  from  a  single 
Cl  measurement.  To  ameliorate  this  problem,  Meier  and  Baladi  estimate  the 
total  unit  weight  and  shear  modulus  of  the  soil  based  on  its  Unified  Soil 
Classification  System  (USCS)  soil  classification,  then  calculate  cohesion  and 
friction  angle  from  two  Cl  measurements  taken  at  different  depths.  This  is 
accomplished  by  simultaneously  solving  two  cavity  expansion  equations  (one  at 
each  depth)  for  the  two  remaining  unknowns. 

4.  The  biggest  drawback  to  the  Meier  and  Baladi  approach  is  that  it 
assumes  the  material  is  homogeneous  between  the  two  measurement  depths.  The 
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methodology  works  fairly  well  for  granular  materials  because  they  are 
reasonably  homogeneous  across  short  distances.  The  properties  calculated  for 
cohesive  materials,  however,  are  not  always  representative  of  the  actual  in 
situ  soil  strength  because  the  in  situ  strength  varies  with  depth  due  to 
changes  in  parameters  such  as  stress  history,  water  content,  and  root  density. 

5.  Farr  (1991)  eliminated  the  need  for  two  Cl  measurements  at  different 
depths  by  estimating  three  of  the  four  unknowns  (shear  modulus,  unit  weight, 
and  friction  angle)  based  on  the  USCS  soil  classification  and  solving  for  the 
remaining  unknown  (cohesion)  using  the  measured  Cl.  The  drawback  to  this 
approach  is  that  it  does  not  make  full  use  of  the  measured  Cl  in  cohesionless 
soils.  Because  the  cohesion  is  not  actually  an  unknown- -it  is  supposed  to  be 
either  identically  zero  or  very  nearly  zero- -there  is  no  need  to  solve  for  it 
using  the  measured  Cl.  The  measured  Cl  is  merely  used  to  confirm  that  the 
three  estimated  soil  properties  (friction  angle,  shear  modulus,  and  unit 
weight)  are  consistent  with  a  soil  having  little  or  no  cohesion.  This 
deficiency  is  addressed  by  the  work  presented  here. 

Purpose 

6.  This  report  describes  the  development  of  a  simple  methodology  for 
estimating  the  unit  weight,  friction  angle,  and  elastic  shear  modulus  of  dry, 
cohesionless  soils  based  solely  on  their  classification  within  the  USCS  and  Cl 
values  obtained  from  the  field.  The  methodology  is  based  on  a  cavity  expan¬ 
sion  formulation  of  the  cone  penetration  problem  and  empirical  correlations 
between  the  mechanical  and  physical  properties  of  soils. 

Scope 

7.  Part  II  presents  an  overview  of  the  theoretical  basis  for  the  pro¬ 
posed  methodology.  Part  III  details  the  algorithms  used  to  estimate  the 
mechanical  properties  of  soils  using  Cl  and  soil  classification.  Part  IV 
presents  a  partial  validation  of  the  methodology  which  is  limited  by  the 
availability  of  suitable  validation  data.  Part  V  contains  a  summary  and  the 
cone lus ions . 
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8.  This  report  also  contains  four  appendixes.  Appendix  A  enumerates 
the  soil  designations  used  in  the  USCS.  Appendix  B  presents  the  mathematical 
derivation  of  the  cavity  expansion  theory  upon  which  the  methodology  is  based. 
Appendix  C  contains  a  listing  of  a  FORTRAN  subroutine  implementing  the  method¬ 
ology.  Appendix  D  contains  a  listing  of  the  notation  usc;d  in  this  report. 


8 


PART  II;  THEORETICAL  FORMULATION  OF  CONE  INDEX 


Background 

9.  There  are  many  correlations  between  cone  penetration  resistance  and 
soil  strength  available  for  use  in  analyzing  foundation  engineering  problems. 
Unfortunately,  very  little  research  has  been  directed  toward  the  more  complex 
problem  of  characterizing  the  surficial  soils  (i.e.,  to  depths  of  1  or  2  ft*) 
that  are  of  interest  in  mobility  applications.  Karafiath  and  Nowatzki  (1978) 
present  a  solution  for  the  penetration  resistance  of  surficial  soils  based  on 
plasticity  theory.  Their  formulation  is  useful  for  investigating  the  effects 
of  cone  geometry  and  penetration  depth,  but  plasticity  theory  does  not  take 
the  stiffness  characteristics  of  the  soil  into  account.  It  is  well  known  that 
soil  compressibility  substantially  affects  penetration  resistance  (e.g.,  Vesic 
1963;  Robeitson  and  Campanella  1983),  so  their  approach  is  not  particularly 
well  suited  to  investigating  the  influence  of  soil  properties  on  penetration 
resistance . 

10.  Rohani  and  Baladi  (1981)  used  the  cavity  expansion  theory  developed 
by  Vesic  (1972)  in  their  formulation  of  cone  penetration  resistance.  Cavity 
expansion  theory  has  been  used  extensively  for  describing  cone  penetration  and 
pile  penetration  resistance  at  depth  (e.g.,  Baligh  1976;  Vesic  1977;  Baldi 

et  al . ,  1981;  Greeuw  et  al.  1988)  because  it  takes  both  the  strength  and  the 
stiffness  of  the  soil  into  account.  Stiffness  is  a  particularly  important 
consideration  in  the  penetration  of  surficial  soils  because  the  soil  is  not 
vertically  constrained  and  will  be  displaced  around  the  cone  and  toward  the 
ground  surface,  thus  lessening  its  apparent  stiffness.  Because  it  takes 
stiffness  into  account,  cavity  expansion  theory  has  been  adopted  as  the  basis 
of  the  methodology  described  here. 


*  A  table  of  factors  for  converting  non-SI  to  SI  (metric)  units  of  measure¬ 
ment  is  presented  on  page  4, 
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Formulation  of  the  Problem 


11.  The  standard  military  cone  penetrometer  consists  of  a  proving  ring, 
a  micrometer  dial,  and  a  30-degree  right  circular  cone  with  a  base  area  of 
0.5  in.^  connected  to  an  aluminum  staff  (see  Figure  1).  While  the  cone  is 
advanced  into  the  soil  by  hand  at  a  more-or-less  constant  rate,  the  micrometer 
dial  is  read,  usually  at  1-in.  penetration  intervals,  to  obtain  the  penetra¬ 
tion  resistance.  The  micrometer  dial  is  calibrated  to  read  the  Cl  directly. 
(The  simplicity  of  the  test  is  self-evident!) 

12.  The  Cl  is  actually  derived  from  the  vertical  force  F^  acting  on 
the  cone  tip  as 


Cl 


n 


(1) 


where  D  is  the  diameter  of  the  cone  tip.  This  equation  converts  the  mea¬ 
sured  force  to  an  index  value  having  units  of  stress. 

13.  The  basic  geometry  of  the  cone  tip  is  shown  in  Figure  2.  The  cone 
tip  can  be  characterized  by  its  base  diameter  D  ,  length  L  ,  and  apex  angle 
2a  (which  are,  of  course,  related  through  L  tan  a  =  D/2  ).  For  the  standard 
WES  mobility  cone,  L  =  1.48  in.,  D  =  0.799  in.,  and  2a  =  30  degrees.  A 
smaller  version  of  the  cone  is  also  available  and  is  often  used  in  particu¬ 
larly  strong  soils  that  would  be  difficult  to  penetrate  by  hand  with  the 
larger  cone.  For  the  smaller  cone  tip,  L  =  0.93  in.  and  D  =  0.5  in.,  thus  the 
base  area  is  0.2  in.^  instead  of  0.5  in.^. 

Stresses  Resisting  Cone  Motion 

14.  If  the  tip  of  the  cone  is  located  at  a  distance  Z-i-L  from  the 
ground  surface,  the  stresses  acting  on  a  finite  frustrum  of  the  cone  at  a 
depth  Z+L-ti  are  shown  in  Figure  2b.  In  that  figure,  a  and  t  are  the 
normal  and  shearing  stresses,  respectively,  that  are  resisting  penetration  of 
the  cone.  Integration  of  these  stresses  over  the  surface  of  the  cone  will 
provide  the  magnitude  of  the  vertical  force  F^  : 
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F, 


(2) 


L 

=  j (a  tana  +  t)  2 nr  dr] 

0 

where  r  —  rj  tan  a  .  The  stresses  a  and  t  are  not  known  a  priori,  so  they 
must  be  determined  from  the  fundamental  properties  of  the  soil. 

15.  It  has  been  observed  in  practice  that  the  cone  tends  to  shear  the 
surrounding  material  during  penetration  and  the  soil  is  observed  to  flow 
around  the  cone  tip  and  towards  the  surface.  This  indicates  that  the  entire 
shearing  strength  of  the  soil  is  being  mobilized,  so  the  normal  and  shearing 
stresses  in  the  soil  can  be  related  through  some  type  of  failure  relationship. 


b  Stresses  on  o  finite  frustrum  of  the  cone 


Figure  2.  The  cone  penetration  problem 
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Assuming  a  Mohr-Coulomb  failure  criterion,  the  normal  and  shearing  stresses 
are  related  as: 


T  =  c  +  o  tan  ((> 


(3) 


where  c  is  the  Mohr-Coulomb  cohesion  intercept  and  (p  is  the  Mohr-Coulomb 
friction  angle. 

16.  To  determine  the  magnitude  of  the  normal  stress  a  acting  on  the 
cone,  the  penetration  of  the  cone  can  be  equated  to  the  expansion  of  a  series 
of  spherical  cavities  and  the  normal  stress  computed  as  the  internal  pressure 
needed  to  maintain  a  slow  expansion  of  those  cavities.  Vesic  (1972)  was  the 
first  to  solve  the  cavity  expansion  problem  for  materials  possessing  both 
cohesion  and  friction.  His  solution  is  summarized  in  Appendix  B. 

17.  For  soils  having  both  cohesive  and  frictional  strength  components 
(i.e.,  c  >  0  and  tan  ^  >  0) ,  the  internal  pressure  needed  to  maintain  a  steady 
expansion  of  a  spherical  cavity  can  be  expressed  as: 


o  =  3  (g  +  c  cot  4i)(  ^  -  c  cot  4) 

\  3  -  sin  4)/  " 


(4) 


(from  Appendix  B)  where  q  is  the  ambient  hydrostatic  stress  existing  prior 
to  the  application  of  the  stresses  at  the  cavity  surface. 


m 


4  sin  4) 
3(1+  sin  4>) 


I 


rr 


Jr 

1  *  A  1^ 


Jr 


G 

c  *  q  tan4> 


and  G  is  the  elastic  shear  modulus. 
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18.  Referring  back  to  the  analogue  of  a  series  of  expanding  spherical 
cavities,  the  hydrostatic  stress  applicable  to  each  of  those  cavities  can  be 
defined  by  q  -  yz  ,  where  7  is  the  total  unit  weight  of  the  soil  and  z  is 
the  depth  below  the  ground  surface.  The  In  situ  stress  corresponding  to  a 
depth  of  Z+h-Ti  (Figure  2a)  is  therefore  given  by 


q  =  {Z  *  L  -  i\)  Y 


(5) 


As  the  stresses  are  integrated  down  the  length  of  the  cone,  this  equation  can 
be  used  to  obtain  the  field  stresses  surrounding  the  corresponding  cavities. 

Calculating  the  Cone  Index 

19.  Combining  Equations  1  through  5  and  completing  the  integration  for 
Fj  ,  the  following  expression  for  Cl  is  obtained: 


CJ  =  6  cl  ^  *  sin<l)U  tang  ♦  tan(l> 
1  3  -  3in4)  /\  tana  ■  tan(J) 


Q 


C  COt(|> 


(6) 


where 


Q  _  [C-t-  (Z-t-Z.)  Y  tan4>]  -  [C->-  ( Z*L)  y  tan(l>  +  (2-/n)  L  y  tan(|>]  (C-*- ZYtan<|>) 

(2-m)  (3-m)  (LYtan4»)^ 


20.  For  the  granular  soils  of  interest  here,  the  cohesion  intercept  is 
zero  and  the  expression  above  simplifies  to 


CJ 


6  C"  ^  3in4)W  tana  tany 
\  3  -  sin(t)/\  tana  •  tan(^ 


(7) 


13 


where 


_  [  (Z-t-L)  Ytan4>]  -  [2Ytan4)  (3-/n)  Ly  tan<|)]  (Zy  tan<t)) 

(2-m)  (3-m)  (Lytan<j>)^ 


21.  Equation  7  expresses  the  Cl  in  terms  of  the  fundamental  soil  prop¬ 
erties  <f>  ,  G  ,  and  7  plus  the  depth  and  geometry  of  the  cone.  Our  goal  of 
relating  Cl  measurements  to  fundamental  soil  properties  has  been  met;  however, 
we  are  still  faced  with  the  problem  of  having  a  single  equation  with  three 
unknowns.  The  problem  can  only  be  made  tractable  by  estimating  two  of  the 
three  unknowns  or  relating  two  of  the  unknowns  to  the  third.  That  is  the  crux 
of  the  methodology  adopted  here  and  the  topic  of  the  next  part. 
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PART  III:  SOIL  PROPERTY  ESTIMATION  FOR  COHESIONLESS  SOILS 


Relating  Shear  Modulus  to  Unit  Weight 

22.  An  estimate  of  the  elastic  shear  modulus  can  be  obtained  using  the 
equations  developed  by  Hardin  and  Richart  (1963).  They  showed  that  for  shear¬ 
ing  strain  amplitudes  less  than  approximately  10"'',  the  elastic  shear  modulus 
in  clean  granular  materials  is  dependent  only  upon  the  void  ratio  of  the  soil 
and  the  ambient  stress  level.  For  materials  with  rounded  grains,  this  rela¬ 
tion  is  expressed  by  the  empirical  equation 

Q  ,  2630(2^17  -e)^  s  (8) 


and  for  soils  with  angular  grains,  it  can  be  expressed  as 

G  =  1230  (2.97  -e)^  .  >0.5 

t  j.  '  O  * 


(9) 


where  G  is  the  shear  modulus  in  psi,  is  the  ambient  stress  level  in 

psi,  and  e  is  the  void  ratio. 

23.  The  void  ratio  can  be  computed  from  the  specific  gravity  and 

dry  unit  weight  7^  of  the  soil  as 


e  = 


Id 


-  1 


(10) 


where  7w  62.4  Ib/ft^  is  the  unit  weight  of  water.  For  most  granular 
soils,  2 . 65  <  Gj  <  2 . 69 ,  so  an  average  value  of  2.67  should  suffice.  If  we 
constrain  the  problem  to  dry  soils,  7d  =  7  «rid  Equation  10  reduces  to 
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e 


166.6 

Y 


(11) 


where  7  is  expressed  in  Ib/ft^. 

Accounting  for  Free  Surface  Effects 

24.  The  cavity  expansion  equations  are  based  upon  the  assumption  that 
the  cavity  expands  in  an  unbounded  medium.  The  process  of  cone  penetration, 
however,  takes  place  in  a  medium  which  is  bounded  above  by  a  free  surface. 

The  existence  of  this  free  boundary  allows  an  upward  flow  of  the  near-surface 
materials  in  the  vicinity  of  the  cone  that  consequently  reduces  the  penetra¬ 
tion  resistance  of  the  cone.  This  free-surface  effect  is  especially  pro¬ 
nounced  in  granular  materials,  but  only  within  about  1  ft  of  the  ground 
surface.  Rohani  and  Baladi  (1981)  present  detailed  discussions  of  the  free 
surface  effect. 

25.  To  account  for  the  free-surface  effect,  Rohani  and  Baladi  postul¬ 
ated  that  the  shear  modulus  G  in  Equations  6  and  7  should  be  replaced  by  an 
apparent  shear  modulus  Gg  that  varies  with  the  penetration  depth  Z  accord¬ 
ing  to  the  equation: 


0.5 


A  + 


1  - 

1  +  Be-P^ 


G 


(12) 


where  A  ,  B  ,  and  0  are  constants  that  reflect  the  cone  geometry  and  must 
be  evaluated  expei imentally .  For  the  standard  WES  cone  operating  in  cohesion¬ 
less  soils,  A  -  0.986,  B  =  100,  and  0  =  0.55  in."^. 

26.  Equation  12  assumes  that  G  is  a  constant.  For  the  surficial 
soils  of  interest  here,  we  can  take  G  to  be  the  shear  modulus  at  the  12-in. 
depth  where  the  free  surface  effects  become  negligible.  At  that  depth,  the 
hydrostatic  stress  is  approximately  1  psi,  and  Equations  8  and  9  conve¬ 

niently  reduce  to 
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(13) 


^  ^  2630(2.17-e)^ 
1  +  e 


and 


Q  ^  1230(2.97  -e)^ 
1  +  e 


(14) 


respectively.  These  equations  uniquely  relate  the  shear  modulus  to  the  dry 
unit  weight  of  the  soil  (by  virtue  of  the  void  ratio  being  uniquely  determined 
by  the  dry  unit  weight) .  If  we  can  also  develop  a  correlation  between  fric¬ 
tion  angle  and  dry  unit  weight,  we  will  have  a  set  of  three  equations  in  three 
unknowns  and  the  problem  will  become  solvable. 

Relating  Friction  Angle  to  Unit  Weight 

27.  Many  researchers  have  developed  correlations  between  friction  angle 
and  dry  unit  weight  for  cohesionless  soils.  The  majority  of  those  correla¬ 
tions  are,  however,  developed  for  just  one  or  two  specific  soils.  Conse¬ 
quently,  they  are  inapplicable  to  the  problem  at  hand.  We  require  a 
correlation  that  applies  to  a  wide  range  of  cohesionless  soil  types,  even  if 
it  has  less  accuracy  than  the  soil-specific  correlations  available  in  the 
literature . 

28.  Figure  3  is  taken  from  the  Naval  Facilities  Engineering  Command 

Design  Manual  7.01  (NAVFAC  1986).  It  depicts  approximate  correlations  between 
friction  angle  and  dry  unit  weight  for  each  of  the  cohesionless  soil  types  in 
the  uses.  This  is  the  missing  piece  of  the  puzzle.  With  this  chart,  we  now 
have  a  tractable  system  of  three  equations  (  Cl  as  a  function  of  G  ,  ^  ,  and 

;  G  as  a  function  of  1  as  a  function  of  7^  )  with  which  we 

can  compute  the  three  unknown  engineering  properties  of  friction  angle,  shear 
modulus,  and  unit  weight  using  only  Cl  measurements  and  the  USCS 
classification. 
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ANGLE  OF  INTERNAL  FRICTION 

VS  density 

(FOR  COARSE  GRAINED  SOILS  ) 


(®^'^MATERIAL 
'  I  TYPE 


^ - 


^  OBTAINED  FROM 
~  EFFECTIVE  STRESS 
FAILURE  ENVELOPES 
APPROXIMATE  CORRELATION 
IS  FDR  COMESWNLESS 
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PLASTIC  FINES 
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DRY  UNIT  WEIGHT  (ypj.PCF 


1.2  i.l  1.0  0.9  Q8  Q73  0.7  065  06  056  05  045  0.4  035  03  025  0.2  0.15 

VOID  RATIO,  e 


0.55  0.5  045  0.4  0.35  0  3  0.25  0.2  0.15 

POROSITY. n 

(0:2.68) 


Figure  3.  Strength  correlations  for  granular  soils 
(Source:  NAVFAC  1986  pg.  7.1-149) 


Estimating  The  Properties  of  a  Given  Soil 


29.  Because  the  three  equations  discussed  above  are  nonlinear,  a  closed 
form  solution  for  the  system  of  equations  does  not  exist.  Instead,  a  numeri¬ 
cal  solution  procedure  will  have  to  be  adopted.  A  simple  trial-and-error 
procedure  can  be  described  as  follows: 

a.  Assume  a  relative  density  and  establish  corresponding  values  of 

and  from  the  appropriate  curve  in  Figure  3. 

b.  Compute  a  void  ratio  e  corresponding  to  using 

Equation  11. 

c.  Compute  a  shear  modulus  G  from  e  using  either  Equation  13 
or  Equation  14  (depending  on  the  angularity  of  the  grains) . 


d.  Using  the  depth  Z  of  the  Cl  measurement,  adjust  G  for  free 
surface  effects  using  Equation  12. 

e.  Substitute  <t>  ,  7d  .  ,  and  Z  into  Equation  7  to  obtain 

the  theoretical  value  of  Cl. 

f.  If  the  measured  Cl  exceeds  the  theoretical  Cl,  choose  a  higher 
relative  density;  otherwise,  choose  a  lower  relative  density. 

£.  Repeat  the  calculations  until  the  theoretical  and  measured  Cl 
values  agree  to  within  a  suitable  tolerance. 

The  values  of  ^  ,  7^  ,  and  G  that  produce  a  theoretical  Cl  equal  to  the 

measured  Cl  are  the  fundamental  soil  properties  being  sought. 

30.  The  procedure  above  should  always  converge  on  a  solution  because 
all  of  the  equations  involved  are  monotonic.  The  number  of  iterations 
required  for  convergence  will,  of  course,  depend  on  the  manner  in  which  the 
iterations  are  performed.  Since  Cl  increases  smoothly  with  increasing  rela¬ 
tive  density  (in  fact,  the  theoretical  Cl  is  very  nearly  a  linear  function  of 
the  logarithm  of  relative  density) ,  a  simple  binary  iteration  approach  will 
suffice . 

31.  Appendix  C  contains  a  FORTRAN  subroutine  implementing  the  procedure 

described  above.  The  subroutine  uses  binary  iteration  to  converge  on  the  Cl 
to  within  one  percent.  To  prevent  convergence  on  totally  unrealistic  relative 
densities,  the  search  for  a  solution  starts  at  a  relative  density  of  -25  per¬ 
cent  and  ends  at  a  relative  density  of  150  percent.  It  would  admittedly  be 
rare  to  find  soils  at  relative  densities  significantly  above  100  percent 
(although  it  would  certainly  be  possible  to  find  soils  at  105  to  110  percent). 
This  wide  range  of  relative  densities  allows  for  soils  that  are  more  dense  in 
nature  than  they  can  be  made  using  standard  laboratory  procedures.  It  also 
provides  some  leeway  for  soils  that  do  not  match  the  NAVFAC  chart  exactly.  If 
a  solution  cannot  be  found  within  the  range  -25  percent  <  <  150  percent, 

the  subroutine  is  exited  via  an  alternate  return. 

32.  Figure  3  has  been  implemented  in  the  subroutine  as  a  tabulation  of 
friction  angle  and  dry  unit  weight  values  for  each  intersection  point  in  the 
chart.  These  values  are  given  in  Tables  1  and  2,  respectively.  Simple  linear 
interpolation  within  these  tables  is  used  to  approximate  the  original  curves. 
In  order  to  have  a  single,  unique  relationship  for  each  soil  type,  it  was 
assumed  a  priori  that  SP  soils  will  fall  to  the  left  of  SM  soils  in  the 
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Table  1.  Friction  angle  versus  relative  density. 


uses 

Relative  Density 

(%) 

Soil 

Type 

0 

25 

50 

75 

100 

ML 

25.9 

28.4 

31.1 

33.4 

36.0 

SP 

26.4 

29.1 

31.9 

34.7 

37.8 

SM 

26.9 

29.7 

32.9 

35.9 

39.1 

SW 

27.1 

30.1 

33.6 

37.0 

40.7 

GP 

27.5 

30.7 

34.6 

38.3 

42.5 

GW 

27.9 

31.4 

36.0 

40.1 

45.0 

Table  2.  Dry  unit  weight  versus  relative  density. 


uses 

Relative  Density 

(%) 

Soil 

Type 

0 

25 

50 

75 

100 

ML 

79.3 

83.9 

88.4 

92.9 

97.9 

SP 

87.4 

92.1 

97.1 

102.3 

107.7 

SM 

95.0 

99.6 

105.0 

110.0 

116.0 

SW 

102,1 

106.4 

112.1 

117.7 

124.3 

GP 

109.0 

114.0 

120.0 

126.6 

133.7 

GW 

117.7 

123.3 

130.6 

138.1 

147.7 
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figure.  Thus,  the  order  of  the  curves  in  Figure  3  is  assumed  to  be:  ML,  SP, 
SM,  SW,  GP,  GW. 

33.  Because  the  shape  of  the  soil  grains  will  seldom  be  known  (in  many 
cases,  the  soil  type  may  not  even  be  known  and  will  have  to  be  assumed),  the 
subroutine  chooses  between  the  two  shear  modulus  equations  based  on  the  com¬ 
puted  void  ratio.  Hardin  and  Richart  (1963)  limit  their  formula  for  rounded 
grains  to  void  ratios  less  than  0.8  and  imply  that  their  formula  for  angular 
grains  only  applies  at  void  ratios  of  0.6  or  greater.  This  is  based  as  much 
on  the  range  of  their  empirical  data  as  anything;  however,  because  soils  with 
rounded  grains  will  not  pack  as  tightly  as  those  with  angular  grains  for  a 
given  relative  density,  it  may  be  assumed  that  soils  with  low  void  ratios  more 
than  likely  nave  rounded  grains,  and  those  with  high  void  ratios  more  than 
likely  have  angular  grains.  This  is  exemplified  by  Figure  4,  which  shows 
minimum  and  maximum  void  ratios  for  two  materials  with  rounded  grains  and  two 
with  angular  grains.  Based  on  these  distinctions,  a  single,  unique  relation¬ 
ship  between  shear  modulus  and  void  ratio  was  obtained  by  using  the  formula 
for  rounded  grains  at  void  ratios  below  0.6,  the  formula  fox  angular  grains  at 
void  ratios  above  0.8,  and  a  weighted  average  of  the  two  for  void  ratios  in 
between. 


A  Caveat  for  Gravel  Soils 


34.  Figure  5  shows  an  idealization  of  a  MCP  penetrating  a  gravelly 
soil.  In  order  for  the  cone  tip  to  advance,  the  individual  gravel  particles 
must  either  fail  in  shear  or  be  pushed  out  of  the  way.  Since  the  MCP  is  hand 
operated,  the  loads  are  insufficient  to  fail  the  gravel  particles;  thus,  par¬ 
ticle  rearrangement  must  control  the  penetration  process.  If  individual 
gravel  particles  are  being  pushed  out  of  the  way  or  rotated  into  more  compact 
orientations,  the  continuum  representation  implied  in  the  cavity  expansion 
theory  is,  strictly  speaking,  invalid.  Great  care  must  therefore  be  used  in 
interpreting  MCPT  results  in  GP  and  GW  soils  using  this  methodology.  The 
further  the  soil  deviates  from  the  continuum  assumptions,  the  greater  the 
potential  error  in  the  soil  property  estimates.  This  means  that  a  well-graded 
gravel  with  a  one -half  in.  maximum  particle  size  is  probably  amenable  to 
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VOID  RATIO 


ROUND  ANGULAR 

COEFFICIENT  OF  ANGULARITY  ,  E 

°  MINIMUM  VOID  RATIO  BASED  ON  MODIFIED 
PROCTOR  COMPACTION  TEST.  EXCEPT  FOR  GLASS  BEADS. 

Figure  4.  Effect  of  particle  shape  on  void  ratio 
(Source;  Holubec  and  D'Appolonia  1972  pg.  308) 


Figure  5,  The  MCPT  in  gravelly  soils 


interpretation  using  this  method.  On  the  other  hand,  a  poorly  graded  gravel 
predominated  by  2-in.  particles  cannot  be  properly  analyzed. 
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PART  IV:  VALIDATION  OF  SOIL  PROPERTY  ALGORITHMS 


Introduction 


35.  Ideally,  validation  of  the  algorithms  presented  in  the  preceding 
chapter  would  be  accomplished  using  MCP,  shear  strength,  and  density  tests 
conducted  in  situ  and  in  close  proximity  to  each  other.  The  latter  require¬ 
ment  is  necessary  to  eliminate,  or  at  least  mitigate,  the  problems  that  can 
arise  from  areal  variability  in  the  soil  properties.  While  in  situ  density 
tests  are  quite  common,  in  situ  strength  tests  are  seldom  performed.  In  the 
dry,  cohesionless  soils  of  interest  here,  an  alternative  would  be  to  substi¬ 
tute  empirical  relations  between  shear  strength  and  dry  density  derived  from 
laboratory  tests  on  remolded  soils  for  the  in  situ  strength  tests.  This  can 
be  done  as  long  as  the  MCPT  and  density  tests  are  performed  in  carefully  com¬ 
pacted  fills  devoid  of  horizontal  layering  and  other  density  variations. 

36.  An  extensive  literature  search  revealed  very  little  data  that  would 
be  suitable  for  validating  the  methodology  presented  in  the  previous  chapter. 
In  the  past,  when  empirical  relationships  based  solely  on  Cl  were  used  for 
almost  all  mobility  modeling,  there  was  no  need  to  directly  measure  fundamen¬ 
tal  soil  strength  properties  such  as  Mohr-Coulomb  friction  angles  and  cohesion 
intercepts.  As  a  result,  such  data  were  rarely  collected  and  seldom  published 
(even  when  they  were  collected).  Thus,  the  search  for  validation  data  is, 
practically  speaking,  constrained  to  experimental  research  programs. 

37.  In  the  research  arena,  there  is  a  plethora  of  data  available  from 
test  programs  conducted  in  calibration  chambers.  Unfortunately,  those  test 
programs  almost  exclusively  simulated  cone  penetration  at  depths  of  many  tr-'s 
of  feet  and  stress  levels  far  in  excess  of  those  assumed  here.  What  little 
research  has  been  conducted  in  surficial  soils  primarily  used  cone  penetration 
resistance  as  the  sole  measure  of  soil  strength,  so  there  are  no  fundamental 
strength  properties  to  compare  against. 

38.  The  authors  managed  to  find  a  limited  amount  of  validation  data  by 
searching  through  WES  reports  and  memoranda  from  the  1960 's  when  a  significant 
amount  of  experimental  research  was  undertaken.  Unfortunately,  all  of  the 
data  was  obtained  in  two  remolded  sands,  both  classified  as  SP  under  the  USCS . 
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With  this  in  mind,  the  remainder  of  this  chapter  presents  a  very  limited  vali¬ 
dation  of  the  soil  property  prediction  methodology. 

Description  of  the  Soils 

39.  Two  sets  of  unpublished  data  on  remolded  sands  were  obtained  from  a 
test  program  conducted  in  the  Large  Blast  Load  Generator  (LBLG)  facility  at 
WES  in  the  mid-1960's.  The  materials  used  in  those  test  programs  were  air-dry 
Cook's  Bayou  and  Reid-Bedford  Model  sands. 

40.  Cook's  Bayou  sand  is  a  uniform,  fine  sand  taken  from  a  borrow  pit 
near  Cook's  Bayou  in  the  vicinity  of  Vicksburg,  Mississippi.  The  material  is 
classified  as  SP  according  to  the  USCS.  Its  minimum  and  maximxim  dry  densities 
are  93.3  and  110.8,  respectively.  Its  specific  gravity  is  2.65,  so  its  maxi¬ 
mum  and  minimum  void  ratios  are  0.77  and  0.49,  respectively.  The  grain  shape 
is  subrounded. 

41.  Reid-Bedford  Model  sand  (hereafter  referred  to  simply  as  Reid- 
Bedford  sand)  is  a  uniform,  fine  sand  obtained  from  a  borrow  pit  near  the  Big 
Black  River  in  the  vicinity  of  Vicksburg,  Mississippi.  The  name  is  derived 
from  the  hydrologic  model  at  WES  in  which  it  was  first  used.  The  material  is 
classified  as  SP  under  the  USCS.  Its  minimum  and  maximum  dry  densities  are 
87.2  and  104.2  pcf,  respectively.  Its  specific  gravity  is  between  2.65  and 
2.66,  so  its  maximum  and  minimum  void  ratios  are  approximately  0.90  and  0.59, 
respectively.  The  grain  shape  is  on  the  border  line  between  subrounded  and 
subangular . 

42.  The  strength  properties  of  both  sands  have  been  documented  in  WES 
reports  by  McNulty  (1965)  and  Kennedy,  Albritton,  and  Walker  (1966).  These 
are  summarized  as  strength-density  curves  in  Figures  6  and  7.  Also  included 
in  those  figures  are  least-square  regression  lines  that  are  plotted  over  the 
range  of  densities  from  minimum  to  maximum  for  each  soil.  The  regression  line 
for  the  Cook's  Bayou  sand  is  given  by 


(t>  =  0.693  -  33.6 


(15) 


and  the  regression  line  for  Reid-Bedford  sand  is  given  by 
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<J)  =  0 . 529  Yd  -  19-3 


(16) 


where  the  dry  unit  weight  is  in  pcf  and  the  friction  angle  <t)  is  in 

degrees . 


Validation  of  the  Methodology 

43.  The  purpose  of  the  LBLG  test  program  from  which  the  validation  data 
were  obtained  was  to  study  the  response  of  buried  arch  structures  to  blast 
loading.  Six  test  specimens,  each  23  ft  in  diameter  and  10  ft  deep,  were 
constructed  to  contain  the  model  arches.  The  specimens  were  built  in  6 -in. 
lifts  using  a  hand- towed  vibrating  plate  to  achieve  a  consistent  specimen 
density.  At  the  completion  of  every  other  lift,  a  series  of  box- density  and 
cone -penetration  tests  were  conducted  at  various  locations  across  the  specimen 
surface.  For  validation  purposes,  tests  conducted  in  close  proximity  to  the 
model  arches  have  been  ignored,  leaving  approximately  50  Cl  profiles  for  each 
specimen.  These  Cl  profiles  consisted  of  readings  taken  at  1-in.  depth 
intervals  for  a  distance  of  12  in.  or  until  the  range  of  the  micrometer  dial 
was  exceeded.  For  our  use,  the  Cl  profiles  have  been  tabulated  at  2 -in.  depth 
intervals  and  an  average  Cl  profile  has  been  obtained  for  each  specimen  by 
averaging  the  cone  indices  at  each  depth.  The  averaged  Cl  profiles  for  the 
six  test  specimens  (designated  here  as  LBLG-1  through  LBLG-6)  are  given  in 
Table  3. 

44.  In  support  of  the  large-scale  testing,  four  small-scale  specimens 
were  built  by  raining  the  sand  into  a  48-in.-diam  mold.  Each  of  these  small 
specimens  was  initially  built  to  a  height  of  18  in. ,  at  which  point  five  cone 
penetration  tests  were  performed.  Following  this,  the  upper  6  in.  of  each 
specimen  was  removed  and  box  density  tests  were  conducted  on  the  top  2  in.  of 
the  remaining  material.  This  provided  dry  unit  weight  data  at  approximately 
the  mid-height  of  the  cone  profiles.  Cone  penetration  data  from  only  one  of 
those  specimens  were  available.  That  data  have  been  averaged  as  well  and  is 
included  in  Table  3  as  specimen  LBLG-7. 

45.  For  each  specimen  in  Table  3,  the  subroutine  listed  in  Appendix  C 
was  used  to  estimate  a  friction  angle  and  a  dry  unit  weight  at  each  tabulated 


27 


Table  3.  Average  cone  profiles  for  the  LBLG  test  specimens. 


Specimen 

Number 

Soil 

Type 

Mean  Dry 
Unit  Weight 
(pcf) 

Cone 

Diameter 
(in. ) 

Cone 
Depth 
(in.  ) 

Mean 

Cone  Index 
(psi) 

LBLG-1 

Cook's  Bayou 

101.4 

0.8 

2 

40 

4 

91 

6 

182 

8 

275 

LBLG-2 

Cook's  Bayou 

107.0 

0.5 

2 

58 

4 

143 

6 

282 

8 

439 

10 

521 

LBLG-3 

Feid-Bedford 

100.5 

0.5 

2 

27 

4 

60 

6 

107 

8 

201 

10 

280 

12 

383 

LBLG-A 

Cook's  Bayou 

107.4 

0.5 

2 

56 

4 

133 

6 

253 

8 

397 

10 

533 

LBLG- 5 

Cook's  Bayou 

101.7 

0.5 

2 

36 

4 

69 

6 

124 

8 

226 

10 

289 

12 

363 

LBLG-6 

Reid-Bedford 

106.5 

0.5 

2 

68 

4 

154 

6 

290 

8 

493 

10 

611 

LBLG -7 

Cook's  Bayou 

107.8 

0.5 

2 

35 

4 

86 

6 

152 

8 

225 

10 

290 

12 

360 
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depth.  The  soil  type  and  Cl  data  input  to  the  subroutine  and  the  friction 
angle  and  dry  unit  weight  estimates  output  from  it  are  presented  in  Table  U. 
For  each  test  specimen,  a  comparison  between  the  soil  properties  measured  in 
the  LBLG  (actually,  the  dry  unit  weights  measured  in  the  LBLG  and  friction 
angles  computed  from  those  data  using  Equations  15  and  16)  and  the  average  of 
the  predicted  soil  properties  are  given  in  Table  5. 

46.  To  better  illustrate  the  amount  of  agreement,  the  predicted  fric¬ 
tion  angles  are  plotted  against  the  measured  friction  angles  in  Figure  8  and 
the  predicted  dry  unit  weights  are  plotted  against  the  measured  dry  unit 
weights  in  Figure  9.  Both  figures  show  very  good  agreement  between  the  pre¬ 
dicted  and  the  measured  values.  There  seems  to  be  a  slight  tendency  towards 
overprediction,  especially  for  the  unit  weights,  but  the  degree  of  overpredic¬ 
tion  is  slight. 

47.  The  data  presented  here  by  no  means  constitute  complete  validation 
of  the  model;  however,  it  does  show  the  potential  predictive  capabilities  of 
the  methodology.  Based  of  these  results,  continued  work,  in  the  form  of  fur¬ 
ther  validation  and  refinement  of  the  model,  is  definitely  warranted. 
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Table  4.  Soil  property  predictions  for  the  LBLG  specimens. 


Soil 

Specimen 

Input  Data 

Output  Data 

Cone 
Width 
(in. ) 

Cone 
Depth 
(in.  ) 

Number 

of 

Tests 

Average 

Cl 

(psi) 

Dry  Unit 
Weight 
(pcf) 

Friction 

Angle 

(deg) 

LBLG-1 

SP 

0.8 

2 

40 

38.7 

4 

91 

■nflH 

37.8 

6 

182 

107.7 

37.8 

8 

44 

275 

106.1 

36.9 

LBLG -2 

SP 

0.5 

2 

58 

42.5 

4 

143 

42.3 

6 

50 

282 

115.1 

42.1 

8 

50 

439 

113.5 

41.1 

10 

44 

521 

110.7 

39.5 

LBLG- 3 

SP 

0.5 

2 

56 

27 

102.3 

34.7 

4 

56 

60 

100.7 

33.9 

6 

56 

107 

99.4 

33.2 

8 

56 

201 

101.3 

34.1 

10 

46 

280 

101.3 

34.1 

12 

46 

383 

102.8 

35.0 

LBLG-4 

SP 

0.5 

2 

50 

56 

HjRHEH 

42.3 

4 

50 

133 

41.6 

6 

50 

253 

113.1 

40.9 

8 

50 

397 

111.7 

40.1 

10 

48 

533 

110.9 

39.7 

LBLG -5 

SP 

0.5 

2 

36 

MWM 

37.6 

4 

69 

35.0 

6 

124 

34.3 

8 

226 

102.8 

35.0 

10 

289 

101.5 

34.3 

12 

363 

102.0 

34.6 

LBLG-6 

SP 

0.5 

2 

26 

68 

■■■ 

■■■ 

4 

26 

154 

6 

26 

290 

B^H 

8 

26 

493 

10 

18 

611 

113.1 

40.9 

LBLG-7 

SP 

0,5 

2 

20 

35 

106.9 

37.3 

4 

20 

86 

106.6 

37.2 

6 

20 

152 

104.5 

35.9 

8 

20 

225 

102.8 

35.0 

10 

20 

290 

101.5 

34.3 

12 

20 

360 

102.0 

34.6 
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Table  5.  Comparison  of  predicted  and  measured  values. 


Soil 

Specimen 

Soil 

Type 

Dry  Unit  Weight  (pcf) 

Friction  Angle  (deg) 

Measured 

Predicted 

Measured 

Predicted 

LBLG-1 

Cook's  Bayou 

101.4 

107.7 

36.7 

37.8 

LB LG-2 

Cook's  Bayou 

107.0 

114.1 

40.6 

41.5 

LBLG-3 

Reid-Bedford 

100.5 

101.3 

33.8 

34.2 

LBLG-4 

Cook's  Bayou 

107.4 

113.1 

40.8 

40.9 

LBLG-5 

Cook's  Bayou 

101.7 

103.0 

34.5 

35.1 

LBLG-6 

Reid-Bedford 

106.5 

115.2 

40.2 

42.1 

LBLG-7 

Cook's  Bayou 

107.8 

104.1 

37.7 

35.7 

Figure  8.  Predicted  friction  angle  versus  measured 

friction  angle 


Predicted  Unit  Weight  (pci) 


Figure  9.  Predicted  dry  unit  weight  versus  measured 

dry  unit  weight 
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PART  V:  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


48.  A  methodology  has  been  presented  for  estimating  the  fundamental 
engineering  properties  of  cohesionless  soils  from  data  collected  with  the  MCP. 
The  theoretical  approach  couples  static  force  equilibrium  with  cavity  expan¬ 
sion  theory  to  model  the  response  of  the  soil  to  penetration  by  a  cone  tip. 

The  mathematical  description  of  the  cone  penetration  process  is  augmented  by 
several  empirical  relationships  that  serve  to  reduce  the  number  of  unknowns 
and  make  the  problem  tractable.  Limited  validation  using  unpublished  data 
collected  in  two  SP  soils  shows  that  the  technique  can  provide  reasonable 
estimates  of  total  unit  weight  and  Mohr-Coulomb  friction  angle.  This  partial 
validation  suggests  that  the  methodology  warrants  further  development. 


Conclusion^  and  Recommendations 


49.  Additional  va' i'  .ion  data  are  needed  before  this  method  can  be 
used  without  reservat?  ^n.  Continued  literature  searches  and  additional  labo¬ 
ratory  and  field  experimentation  are  needed  to  further  validate  and  refine  the 
methodology.  In  order  to  be  usable,  validation  data  must  include  in  situ  mea¬ 
surements  of  unit  weight,  shear  strength,  shear  stiffness,  and  cone  penetra¬ 
tion  resistance- -all  taken  in  close  proximity  to  one  another. 

50.  Future  research  should  be  directed  towards  refining  the  existing 
methodology  for  cohesionless  soils  and  extending  the  theory  to  cohesive  soils. 
It  is  hoped  that,  eventually,  a  system  can  be  devised  that  is  applicable  in 
all  -.oil  types.  One  approach  would  be  to  collect  several  methodologies  such 
as  the  one  described  here  and  embed  them  within  a  knowledge -based- system 
(KBS).  A  KBS  (which  includes,  but  is  not  limited  to,  so-called  "expert  sys- 
•“ams")  is  a  computer  program  that  combines  knowledge  about  a  specific  domain 
ind  reasoning  abilities  to  arrive  at  conclusions.  In  the  domain  of  soil- 

property  estimation,  a  KBS  might  combine  cone  penetration  and  soil  description 
data  supplied  by  the  user  with  a  collection  of  theoretical  and  empirical 
models  to  arrive  at  estimates  of  the  soil  properties.  The  KBS  would  choose 
the  appropriate  model (s) ,  infer  or  ask  for  any  missing  data,  and  sort  out  the 
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often  conflicting  answers  that  the  models  produce.  The  research  described 
here  is  a  first  step  towards  producing  such  a  system. 
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APPENDIX  A:  UNIFIED  SOIL  CLASSIFICATION  SYSTEM 


Table  A1 

Detailed  description  of  the  USCS  soil  classes 


uses  Soil 
Classification 

GW 

GP 

GM 

GC 

SW 

SP 

SM 

SC 

ML 

CL 

OL 

MH 

CH 

OH 

PT 


_ Description _ 

Well  graded  gravels,  gravel-sand  mixtures,  little  or 
no  fines,  wide  range  of  grain  sizes 

Poorly  graded  gravels,  gravel -sand  mixtures,  little  or 
no  fines,  narrow  range  of  grain  sizes 

Silty  gravels,  poorly  graded  gravel-sand-silt  mix¬ 
tures,  nonplastic  fines  or  fines  with  little 
plasticity 

Clayey  gravels,  poorly  graded  gravel -sand- clay  mix¬ 
tures,  plastic  fines 

Well  graded  sands,  gravelly  sands,  little  or  no  fines, 
wide  range  of  grain  sizes 

Poorly  graded  sands,  gravelly  sands,  little  or  no 
fines,  narrow  range  of  grain  sizes 

Silty  sands,  poorly  graded  sand-silt  mixtures,  non¬ 
plastic  fines  or  fines  with  low  plasticity 

Clayey  sands,  poorly  graded  sand-clay  mixtures,  plas¬ 
tic  fines 

Inorganic  silts  and  very  fine  sands,  rock  flour,  silty 
or  clayey  fine  sands  with  slight  plasticity 

Inorganic  clays  of  low  to  medium  plasticity,  gravelly 
clays,  sandy  clays,  silty  clays,  lean  clays 

Organic  silts  and  organic  silty-clays  of  low 
plasticity 

Inorganic  silts,  micaceous  or  diatomaceous  fine  sandy 
or  silty  soils,  elastic  silts 

Inorganic  clays  of  high  plasticity,  fat  clays 

Organic  clays  of  medium  to  high  plasticity 

Peat  and  other  highly  organic  soils 
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APPENDIX  B:  DERIVATION  OF  THE  CAVITY  EXPANSION  EQUATIONS 


1.  The  problem  of  the  quasi -static  expansion  of  a  spherical  cavity  in 
an  unbounded,  compressible,  elastic-plastic  medium  obeying  the  Mohr-Coulomb 
failure  criterion  was  first  solved  by  Vesic  (1972).*  His  solution  is  summa¬ 
rized  here  for  the  convenience  of  the  reader. 

2.  The  objective  is  to  determine  the  radial  stress  Op  at  the  cavity 
surface  necessary  to  maintain  a  slow  expansion  of  the  cavity.  Referring  to 
Figure  B1 ,  stress  equilibrirun  for  an  element  located  at  some  distance  r  from 
the  center  of  the  cavity  can  be  written  as 


(Bl) 


where  is  the  radial  stress  and  ag  is  the  circumferential  stress.  If  it 

is  assumed  that  the  material  immediately  beyond  the  cavity  is  in  a  state  of 
incipient  failure  (i.e.,  limit  equilibrium  is  invoked),  then  the  radial  and 
circumferential  stresses  can  be  related  through  the  failure  criterion  as 

-  Og  =  (o^  +  08)sin<J>  +  2  c  cos4>  (B2) 


where  c  is  the  Mohr-Coulomb  cohesion  intercept  and  (f>  is  the  Mohr-Coulomb 
friction  angle.  Combining  Equations  Bl  and  B2 ,  the  governing  differential 
equation  for  the  radial  stresses  around  the  cavity  becomes 


dz 


^  sin  ^  *  2  c  cos  4> 
r  \  1  +  sin  4) 


(B3) 


This  equation  must  be  solved  subject  to  the  boundary  conditions  shown  in  the 
figure . 

3.  Equation  B3  is  a  first-order  differential  equation  of  the  form 


*  References  cited  in  this  appendix  are  listed  in  the  References  section  at 
the  end  of  the  main  text. 
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where  Cj  is  an  integration  constant.  Applying  the  boundary  condition  a^.  - 
Pu  at  r  “  and  completing  the  integration: 


a 


r 


(p.. 


C  cot  4>) 


4aln  A 
i  ♦  sin^ 


c  cot  <j> 


(B6) 


This  expression  relates  the  radial  stresses  in  the  vicinity  of  the  cavity  to 
the  instantaneous  size  of  the  cavity  and  the  internal  cavity  pressure. 

4.  To  determine  the  pressure  p„  needed  to  expand  the  cavity  to  some 
radius  ,  Vesic  first  applied  the  continuity  condition  that  the  change  in 

volume  of  the  cavity  must  be  equal  to  the  change  in  volume  of  the  elastic  zone 
plus  the  change  in  volume  of  the  plastic  zone: 

rI  =  rI  -  (Rp-Up)^  *  {R^-R^)A  (B7) 


where  Up  is  the  radial  displacement  of  the  boundary  betwee '  the  elastic  and 
plastic  zones,  A  is  the  average  volumetric  strain  in  the  plastic  zone,  and 
it  is  assumed  that  the  initial  cavity  radius  is  zero.  From  elastic  theory, 
the  displacement  of  the  elastic/plastic  boundary  is  proportional  to  the 
increase  in  radial  stress  above  the  initial  ambient  pressure  q  : 


u 


p 


(B8) 


Substituting  Equation  B8  into  Equation  B7  and  ignoring  higher  powers  of  Up  : 


1  ♦  A 


Pv 


(B9) 


Substituting  r  -  i?p  in  Equation  B6  provides  an  expression  for  : 
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(BIO) 


o 


p 


(Py  +  c  cot  4>) 


Up] 


C  cot  <j) 


Stress  equilibrium  requires  that 


(p„  +  c  cot  <t») 


4 

1  ♦  41114 


c  cot  (|> 


(BID 


Substituting  Equations  BIO  and  Bll  into  Equation  B9 ,  the  following  expression 
for  the  ratio  of  the  plastic  zone  radius  to  Che  cavity  radius  is  obtained: 


vr 


-T-J-  -  //c  +  gtan<J)\ 

j  3  cos  <J>  \ 

V  \  G  ) 

\  3  -  sin  (J)/ 

+  A 


(B12) 


Noting  that  for  A  small 


Vl  ♦  A  =  1 


and  for  <f>  <  45° 


3  cos  4> 
3  -  sin 


Equation  B12  can  be  simplified  and  rearranged  to  produce  the  approximation 


Ru  N  1  *  I,A 


(B13) 


in  which  1^  ,  the  rigidity  index,  is  equal  to  the  ratio  of  the  shear  modulus 
of  the  soil  to  its  initial  strength; 
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c  *  q  tan  4> 


5.  If  the  material  is  assumed  to  be  incompressible  (i.e.,  A  -  0) ,  Equa¬ 
tion  B13  reduces  to 


Ru 


=  {: 


q  tan 


4> 


(B14) 


Substituting  this  equation  back  into  Equation  BIO,  we  arrive  at  the  required 
expression  for  the  stress  at  the  cavity  surface  (i.e.,  at  r  —  ); 


o  =  3  (  (7  +  c  cot4i  )  1 1  *  I  -  lit*)  .  ^  cot*  (B15) 

\  3  -  sin4>/ 

This  is  used  in  Chapter  2  to  derive  an  expression  for  the  normal  stresses 
acting  on  the  cone  surface. 
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APPENDIX  C: 


A  FORTRAN  IMPLEMENTATION  OF  THE  METHODOLOGY 


c  —  I - - - 

C  Binary  Iteration  Routine  for  Solving  Cone  Index  Equation 

C 

C  INPUT  v/ARIABLES 

C 

C  Width  -  Maximum  cone  diameter  in  inches  (real) 

C  uses  -  uses  soil  classification  (2-character  string) 

C  Z  -  Cone  index  measurement  depth  in  inches  (real) 

C  Cl  -  Measured  cone  index  in  psi  (real) 

C 

C  OUTPUT  VARIABLES 

C 

C  Density  -  Relative  density  in  percent  (real) 

C  Gamma  -  Dry  unit  weight  in  pcf  (real) 

C  Phi  -  Friction  angle  in  degrees  (real) 

C  G  -  Elastic  shear  modulus  in  psi  (real) 

C  *  -  Alternate  return  for  error  condition 

C 

C  LOCAL  VARIABLES 

C 

C  Soil  -  Soil  type  index  into  Angles  and  Gammas  arrays  (integer) 

C  Lower  -  Density  index  into  Angles  and  Gammas  arrays  (integer) 

C  Angles  -  Array  of  friction  angles  for  each  soil  type  (real) 

C  Gammas  -  Array  of  unit  weights  for  each  soil  type  (real) 

C  Symbol  -  Array  of  USCS  soil  classification  symbols  (character*2) 

C  E  -  Void  ratio  (real) 

C  AngleG  -  Shear  modulus  for  angular  soils  in  psi  (real) 

C  RoundG  -  Shear  modulus  for  rounded  soils  in  psi  (real) 

C  Predict  -  Computed  cone  index  in  psi  (real) 

C 

C  . 

Subroutine  SOLVER  (USCS , Z ,CI .Density .Gamma , Phi , G ,*) 

Integer  Soil. Lower 

Real  Angles( 5 , 6) ,Gammas( 5 . 6) 

Real  Width , Z , Cl . Phi . G , Gamma . E . Density 

Real  Increment , Ratio . Delta , Fuzz .AngleG . RoundG , Predict 

Character  USCS*2 , Symbol*2 (2 . 6) 


Data  Symbol 

/  'gw' 

'gP' 

'sw' 

.  'GW' 

% 

1 

.  'GP' 

t 

2 

,  'SW' 

» 

3 

'  sm' 

,  'SM' 

4 

'  sp ' 

,  'SP' 

t 

5 

'  ml ' 

,  'ML' 

/ 

Data  Angles 

/  27.9 

,  31.4 

,  36.0 

,  40.1 

,  45.0 

1 

27.5 

,  30.7 

,  34.6 

,  38.3 

,42.5 

2 

27.1 

,  30.1 

,  33.6 

,  37.0 

,  40.7 

3 

26.9 

,  29.7 

,  32.9 

,  35.9 

,  39.1 

4 

26.4 

,29.1 

,  31.9 

,  34.7 

.  37.8 

5 

25.9 

,  28.4 

,  31.1 

,  33.4 

,  36.0 

/ 

Data  Gammas 

/117.7 

,123.3 

,130.6 

,138.1 

,147.7 

1 

109.0 

,114.0 

,120.0 

,126.6 

,133.7 

» 
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CM  ro  m 


102.1  .106.4  ,112.1  .117.7  .124.3  . 

95.0  ,  99.6  .105.0  .110.0  .116.0  , 

87.4  ,  92.1  .  97.1  .102.3  .107.7  , 

79.3  ,  83.9  .  88.4  .  92.9  .  97.9  / 

C  ...  convert  soil  classification  symbol  to  an  array  index 

Do  1  J-1.6 
Do  1  1-1,2 

1  If  (uses. EQ. Symbol ( I, J))  Soil  -  J 


C  ...  initialize  relative  density  and  zero  tolerance  for  convergence 


Density  —  -25.0 
Increment  =  25.0 
Fuzz  -  0.01*CI 


C  ...  convert  relative  density  guess  to  an  array  index 
2  Lower  -  MAX(  0  ,  MIN  (  INT(Density/25)  ,  3  )) 

C  ...  use  linear  interpolation  to  compute  phi  and  gamma  from  density 
Ratio  -  (Density-25*Lower)/25 

Delta  -  Angles(Lower+2 ,Soil)  -  Angles(Lower+l , Soil) 

Phi  -  Angles (Lower+1 , Soil)  +  Delta*Ratio 

Delta  -  Gammas(Lower+2 ,Soil)  -  Gammas(Lower+l , Soil) 

Gamma  -  Gammas (Lower+1 , Soil)  +  Delta*Ratio 


C  ...  compute  void  ratio  from  gamma  and  shear  modulus  from  void  ratio 

E  -  (2.68  *  62.4  /  Gamma)  -  1.0 

If  (E.LE.0.6)  Then 

G  -  2630  *  (2.17-E)**2  /  (1.0+E) 

Else  If  (E.GE.0.8)  Then 

G  -  1230  *  (2.97-E)**2  /  (1.0+E) 

Else 

RoundG  -  2630  *  (2.17-E)**2  /  (1.0+E) 

AngleG  -  1230  *  (2.97-E)**2  /  (1.0+E) 

G  -  (E-0.6)/0.2*RoundG  +  (0 . 8-E)/0 . 2*AngleG 
End  If 

C  ...  compute  cone  index  corresponding  to  phi,  g,  gamma,  and  depth 
Predict  -  C0NE(Phi , G , Gamma , Z .Width) 

C  ...  compare  to  measured  cone  index  and  adjust  density  guess  as  needed 

If  (ABS(Predict-CI) .LT.Fuzz)  RETURN 
If  (Predict. GT. Cl)  Then 

If  (Density . LE .- 25 . 0)  Then 
RETURN  1 
Else 

Increment  -  0 . 5*Increment 
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Density  —  Density  -  Increment 
End  If 
Else 

If  (Density . GE . 150 . 0)  Then 
RETURN  1 
Else 

Density  =  Density  +  Increment 
End  If 
End  If 

C  ...  compute  another  cone  index  estimate  using  the  new  relative  density 
Go  To  2 
Ei»d 
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c 

c 

c 

c 

c 

c 

c 
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Mathematical  Model  for  Computing  Cone  Index 


Phi 

G 

Gamma 

Z 

Width 


INPUT 

Friction  angle  in  degrees  (real) 

Elastic  shear  modulus  in  psi  (real) 

Dry  unit  weight  in  pcf  (real) 

Cone  index  measurement  depth  in  inches  (real) 
Maximum  cone  diameter  in  inches  (real) 


Function  CONE (Phi . G , Gamma , Z .Width) 

Parameter  (Alpha=15 . 0 , A=0 . 986 , B=100 . 0 , Beta=0 . 55 ) 

Double  Precision  GammaTanPhi .Terml ,Term2 .Omega 
Real  Phi , G . Ganuna . Z , Pi , TanAlpha . SinPhi .TanPhi , H , M , Ga 

C  ...  compute  cone  height  from  width  and  apex  angle 

Pi  -  ATAN2(0.0, -1.0) 

TanAlpha  =  TAN(Pi*Alpha/180) 

H  -  0 . 5*Width/TanAlpha 

C  ...  compute  trig  functions  involving  friction  angle 

SinPhi  -  SIN(Pi*Phi/180) 

TanPhi  -  TAN(Pi*Phi/180) 

GammaTanPhi  -  (Gamma/17 28 )*TanPhi 

C  ...  compute  the  exponent  of  the  shear  modulus 

M  -  4./3.  *  SinPhi/(l.+SinPhi) 

C  ...  compute  the  omega  term 

Terml  >=  ( (Z+H)*GammaTanPhi)**(3-M) 

Term2  =  (Z+(3-M)*H)  *  GammaTanPhi  *  (Z*GammaTanPhi)**(2-M) 

Omega  =  (Terml-Term2)  /  (2-M)  /  (3-M)  /  (H*GammaTanPhi)**2 

C  ...  adjust  the  shear  modulus  for  free  surface  effects 

Ga  -  0.5  *  G  *  (A  +  ( 1 . -B*EXP( -Beta*Z) )  /  ( 1 . +B*EXP( -Beta*Z) )  ) 

C  ...  compute  and  return  the  cone  index 


CONE  -  6  *  Ga**M  *  Omega 

*  (l.+SinPhi)  /  (3. 

*  (TanAlpha+TanPhi) 

-SinPhi) 

/  (TanAlpha*TanPhi ) 

Return 

End 
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APPENDIX  D: 


NOTATION 


A 

fi 

c 

Cl 

Cl 

Dr 

E 

e 

Fz 

G 

Ga 

Gs 

Ir 

L 

m 

n 

Pu 

q 

r 

Ru 

Rp 

z 

a 

P 

A 

7 

7d 

»7 

a 

Oa 


Empirical  constant  used  to  calculate  apparent  shear  modulus 
Empirical  constant  used  to  calculate  apparent  shear  modulus 
Mohr-Coulomb  cohesion  intercept 
Constant  of  integration 

Cone  index  (cone  penetration  resistance) 

Relative  density  of  the  soil 

Coefficient  of  angularity  of  the  soil  grains 

Void  ratio  of  the  soil 

Measured  vertical  force  on  the  cone 

Initial  elastic  shear  modulus  of  the  soil 

Apparent  shear  modulus  (adjusted  for  free-surface  effects) 
Specific  gravity  of  the  soil 
Rigidity  index  of  the  soil 
Length  of  the  cone  tip 

Shear  modulus  exponent  in  the  cone  index  equation 
Porosity  of  the  soil 

Pressure  inside  a  spherical  cavity  of  arbitrary  radius 
Overburden  stress  at  an  arbitrary  depth 
Arbitrary  radius  of  the  cone  tip 

Arbitrary  radius  of  an  expanding  spherical  cavity 
Arbitrary  radius  of  the  plastic  zone  around  the  cavity 
Arbitrary  displacement  of  the  elastic-plastic  boundary 
Depth  of  cone  penetration 
Apex  angle  of  the  cone  tip 

Empirical  constant  used  to  calculate  apparent  shear  modulus 

Volumetric  strain  in  the  vicinity  of  the  cavity 

Total  unit  weight  of  the  soil 

Dry  unit  weight  of  the  soil 

Unit  weight  of  water 

Arbitrary  vertical  distance  from  the  cone  tip 
Arbitrary  normal  stress  acting  on  the  cone  tip 
In  situ  confining  pressure 

Pressure  inside  an  expanding  spherical  cavity 
Arbitrary  radial  stress 
Arbitrary  circumferential  stress 
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r 


Arbitrary  shear  stress  acting  on  the  cone  tip 
Mohr-Coulomb  friction  angle 

Term  used  to  simplify  the  cone  index  equation 


4> 
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